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1. Introduction 


The mass composition of ultra-high energy cosmic rays (UHECR) inducing extensive air 
showers of energies above 10*^ eV can be measured by fluorescence detectors on average ba¬ 
sis. The measurement of the depth of shower maximum (Xmax) is compared with Monte Carlo 
(MC) predictions in such cases i@i- The fluorescence technique provides a precise measurement 
of Xmax> but a large systematic uncertainty remains in the determination of the mass composition 
of UHECR. The obstacle comes predominantly from different predictions of hadronic interaction 
models that are extrapolated from accelerator energies to energies larger by a few orders of magni¬ 
tude in the center of mass system. The mass composition of UHECR remains uncertain, and even 
unknown at the highest energies where a steep decrease of the flux is observed [^,^. The very low 
statistics of events collected by fluorescence detectors is due to their low duty cycle. 

Assuming a small number of primaries to be present in UHECR, the most probable fractions 
of these primaries were inferred in |Q]. The measured distributions of X^ax were compared with 
Xmax distributions by combining MC distributions of the assumed primaries. Earge differences in 
the results were found among the hadronic interaction models. Also, a degeneracy of solutions with 
similar probability can be expected as, generally, there are many combinations of MC distributions 
of the individual primaries which describe the measured distributions similarly well. 

Whereas the fluorescence technique measures the longitudinal profile of the electromagnetic 
component of the shower, a measurement of the number of muons on the ground (N^) can pro¬ 
vide an independent way to infer the mass composition of UHECR. Muon detectors have a 100% 
duty cycle and, when a good resolution of is achieved (<10%), even better separability of the 
individual primaries can be achieved than in the case of the fluorescence technique. However, 
there is a lack of in MC simulations when compared with the measured data i 0, i- The 
underestimation of muon production is usually characterized in terms of a muon rescaling fac¬ 
tor. Moreover, a stronger relationship between and the shower energy than between X^ax and 
the shower energy [^] makes the situation more difficult and an independent measurement of the 
shower energy is needed for composition studies using N^. Therefore a simple comparison of the 
measured distributions of muons with MC predictions would be complicated. 

A combined measurement of UHECR showers with the fluorescence technique (X^ax and 
shower energy) and muon detectors (N^) could be a more sucessfull way to determine the mass 
composition of UHECR. In the previous studies, the detected muon and electromagnetic signals 


were utilized to determine the average mass number of a set of air showers, see e.g. [10]. There are 
also methods estimating the spread of masses in the UHECR primary beams via correlation of 
and Xmax [Oil frorn signals in muon and electromagnetic detectors [l^j. However, it needs to 
be mentioned that the two cuiTently operating experiments do not yet directly measure the muonic 
component for showers with zenith angles below 60°. 

Here we present another method to determine the fractions of the assumed primaries in which 
the rescaling of (R^) can be achieved simultaneously with a single fit of a combined measure¬ 
ment of Xmax and N^. Eor this purpose we use MC showers generated with two hadronic interaction 
models tuned to the EHC data (run I-II). 

In the next section, the generated showers that were used in this work are described. The 
method to determine the fractions of the assumed primaries and is introduced in Section 
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Section ^ contains applications of the method on several examples, testing its performance. The 
work is summarized in the last section. 


2. Simulated Showers 


For this work we have simulated ~ 10^ showers with the CONEX 4.37 generator |]3, 14] for 
p and Fe primaries with fixed energy 10^^’^ eV and for each of the two hadronic interaction models 
(QGSJet 11-04 [15] and EPOS-LHC [|^, [^). The zenith angles (0) of showers were distributed 
uniformly in cos^ 0 for 0 in < 0, 60° >. 

Muons with threshold energy 300 MeV at ~1400 m a.s.l. were included to calculate N^. 
Electromagnetic particles of energies above 1 MeV formed the longitudinal profile (dependence of 
fhe deposifed energy on fhe afmospheric depfh), from which Xmax was fiffed wifh fhe Gaisser-Hillas 
function by fhe program CONEX. 

For each shower we used Gaussian smearing of X^ax and wifh a variance equal fo a(Xniax) 
and a(V^), respectively. These smearings imifafe fhe defector resolutions. We adopfed a correcfion 
for fhe aftenuafion of wifh zenifh angle due fo fhe differenf amounf of afmosphere penefrafed 
by fhe air shower before if reaches fhe ground. The correcfion was made using a polynomial of 3’^'* 
order in cos^0 for each model of hadronic interacfions. An equally mixed composition of p and 
Fe was considered for fhis purpose. 




Figure 1: Disfribufions of X^ax paramefrized wifh Gumbel funcfions g,. Showers inifiafed wifh p 
(blue) and Fe (red) primaries were simulafed wifh QGSJef 11-04 (leff) and EPOS-FHC (righf) for 
CT(Xmax) = 20 g/cm^. 


3. Method 

For bofh hadronic inferacfion models, we pai'amefrized X^ax disfribufions wifh Gumbel func¬ 
tions gi (see Fig. for bofh primaries i = p, Fe. We also paramefrized fhe dependence of mean 
>, on Xmax wifh quadrafic funcfions in X^ax denofed as < > (see Fig. again for 

bofh primaries. We infroduced fhe rescaling factor of < > fo incorporafe info fhe mefhod 
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Figure 2: Dependence of mean on parametrized with quadratic functions <N‘^>. Showers 
initiated with p (blue) and Fe (red) primaries were simulated with QGSJet 11-04 (left) and EPOS- 
LHC (right) for a(Xmax) = 20 g/cm^ and o{N^)/N^ = 10%. Only bins of X^ax with more than 30 
showers were considered in the quadratic fits. 


the case when a rescaling of < > obtained from MC is needed to fit the measured < >. 

Then, for a combination of two primaries with fractions fi,Y.fi= > is given as 

<Nyi>=Y.[wi<Nl>)R^ (3.1) 


where the weights w, are expressed as 


W; = 


fi ■ gi 


Lifj-gjY 


(3.2) 


For each bin of Xmax^ < >=< > (Xmax) is calculated as the weighted average of 

< N‘^ >=< N'^ > (Xmax) rescaling both <N'^> with the same factor R^. The weights w,- = w,(Xniax) 
reflect the relative contribution of each individual primary with relative fraction ft in each bin of 
Xmax according to gi = g;(Xmax)- 

Thus, any given dependence of < > on X^ax; which is similar to the dependence of com¬ 

bined proton and iron showers, can be fitted with the two-parameter (/p and R^) fit. The Fe fraction 
is obtained afterwards as /pe = 1 — /p- 

An example of application of the method to the mixed composition of showers initiated with 
50% p and 50% Fe is shown in Fig. The fitted dependence of < > on X^ax (black points) 

is shown with the gray dashed line. The hadronic interaction model EPOS-LHC was considered. 
Starting from the lowest values of Xmax> < Nji > matches < X^® > to about 650 g/cm^ where a 
transition towards < X^ > begins. It continues up to about 800 g/cm^ where < X^^® > starts to 
match < X^ >. 


4. Application of Method 

In this section we show basic examples of the present method and rough estimates how accu¬ 
rately the primary fractions of p and Pe and the muon rescaling factor can be determined. In the 
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Figure 3: Example of a fit (gray dashed line) for a set of showers composed of 50% p and 
50% Fe (black points). Showers were generated with EPOS-LHC for a(Xmax) = 20 g/cm^ and 
(7{N^)/N^ = 10%. The individual p and Fe showers are shown with light blue and light red points, 
respectively. < > is depicted with blue open points and < > with red open points. The tit 

with gi and < N‘^> parametrized for EPOS-EHC is shown with the gray dashed line. 


following, we assumed the detector resolutions to be a(Xmax) = 20 g/cm^ and a{N^)/N^ = 10%. 
We considered 11 combinations of mixed compositions of p and Fe with fractions in steps of 10% 
for both hadronic interaction models. Each of these compositions was reconstructed with each of 
the two parametrizations obtained for the two hadronic interaction models. Additionally, we re¬ 
constructed an example dataset with parametrizations of each of the two models to have another 
assessment of the method with respect to the two models. 

On the left panel of Fig. ^ a comparison of the fitted and the true proton fraction is shown 
for QGSJet 11-04 (blue) and EPOS-EHC (red). The reconstructed /p of showers generated with a 
different model than that used for the parametrization of gi and < > are depicted by open black 

markers. Scenario 1 (2) corresponds to showers produced with QGSJet 11-04 (EPOS-EHC) and 
fitted with the parametrizations obtained from EPOS-EHC (QGSJet 11-04) showers. 

When the same hadronic interaction model is used for the generation of showers and parametriza¬ 
tion of gi and <N‘^>, the proton fraction is reconstructed within a few % of the true value. How¬ 
ever, in cases of Scenario 1 and 2, the difference between the fitted and the true proton fraction 
increases with the spread of primary masses of the selected composition. It reaches values up to 
~ 20 %. 

On the right panel of Fig. Q, the muon rescaling is plotted for different true proton fractions. 
The rescaling factor is found to be within a few % to 1 (precision of the method), when the same 
models were used for parametrization and generation (red and blue). Black points correspond to 
the relative difference of for showers generated with QGSJet 11-04 and EPOS-EHC, which is 
about 6%. 
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Figure 4: Left panel: Comparison of fitted (/p) and true proton fraction (True /p). Right panel: 
Rescaling of depending on the true proton fraction. Perfect matches between fitted and true 
values are shown with dotted gray lines. Smearing with a(Xmax) = 20 g/cm^ and o{N^)/N^ = 10% 
was considered. Scenario 1 (2) corresponds to QGSJet 11-04 (EPOS-LHC) showers fitted with 
parametrizations of g, and < N‘^> from EPOS-LHC (QGSJet 11-04). 
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Eigure 5: An example of data (black points) fitted with QGSJet 11-04 (blue) and EPOS-EHC (red) 
parametrizations. 


As another check, we created example data from 5000 p and 5000 Ee showers produced with 
QGSJet 11-04. Eor each shower, we scaled by a factor 1.3 and increased A^ax by 7 g/cm^. Note 
that EPOS-EHC generates showers with deeper A^ax than QGSJet 11-04 by about 14 g/cm^ on 
average. These example data were fitted with parametrizations of both models (see Eig. Both 
fits describe the example data similarly well giving different proton fractions and muon rescaling 
factors that are shown in Tab. |i|. The difference of fp is about 15% when different parameterizations 
(Eigs. [T],^ based on two most recent models are used. The ratio of for the two models reflects 
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again that EPOS-LHC produces about 6% more muons than QGSJet 11-04 on average. 


Table 1: Fitted parameters for example data. 


Model 

/p [%] 

Rh 

QCSJet 11-04 

41 ±2 

1.297 ±0.004 

EPOS-FHC 

56 ±2 

1.216 ±0.004 


5. Conclusions 

A method of simultaneously obtaining the primary fractions and the muon rescaling factor 
from Xmax and of UHECR was presented. Simulated showers with two models of hadronic 
interactions tuned to EHC data (run I-II) were used. The precision of the method was tested with 
different combinations of p and Fe primaries and with example data. The primary fractions and 
the muon rescaling factor can be determined within a few %. The difference of the proton fraction 
reconstructed with the two parameterizations based on the two models of hadronic interactions was 
observed below 20%. The muon rescaling factor reflected the relative difference (around 6%) in 
the average muon shower size of the two models of hadronic interactions. 

Acknowledgements 

This work is funded by the Czech Science Foundation grant 14-17501S. 

References 

[1] A. Aab et al. (The Pierre Auger Collaboration), Phys. Rev. D 90 (2014) 122005. 

[2] R. U. Abbasi et al. (Telescope Array Collaboration), Asfroparf. Phys. 64 (2014) 49-62. 

[3] J. Abraham et al. (The Pierre Auger Collaboration), Phys. Rev. Lett. 101 (2008) 061101. 

[4] T. Abu-Zayyad et al. (Telescope Array Collaboration), A/?/ 768 (2013) LI. 

[5] A. Aab et al. (The Pierre Auger Collaboration), Phys. Rev. D 90 (2014) 122006. 

[6] A. Aab et al. (The Pierre Auger Collaboration), JCAP 8 (2014) 019. 

[7] B. Kegl for the Pierre Auger Collaboration, Proc. of the 33rd ICRC 2013, arXiv: 1307.5059 
[astro-ph.HE]. 

[8] G. R. Farrar for the Pierre Auger Collaboration, Proc. of the 33rd ICRC 2013, arXiv:1307.5059 
[astro-ph.HE]. 

[9] J. Mathews, Astropart. Phys. 22 (2005) 387-397. 

[10] P. Luczak et al., Proc. of the 33rd ICRC 2013, arXiv: 1308.2059 [astro-ph.HE]. 

[11] P. Younk & M. Risse, Astropart. Phys. 35 (2012) 807-812. 

[12] J. Vicha et al., Astropart. Phys. 69 (2015) 11-17. 


7 









Combined Analysis ofN^ and X^^y^for UHECR showers 


Jakub Vicha 


[13] T. Bergmann et al., Astropart. Phys. 26 (2007) 420-432. 

[14] T. Pierog et al., Nucl. Phys. B. - Proc. Siippl. 151 (2006) 159-162. 

[15] S. S. Ostapchenko, Phys. Rev. D 83 (2011) 014018. 

[16] K. Werner, F. M. Liu and T. Pierog, Phys. Rev. C 74 (2006) 044902. 

[17] T. Pierog and K. Werner, Nucl. Phys. B - Proc. Suppl. 196 (2009) 102-105. 

[18] E. J. Gumbel, Statistics of extremes, Dover Publications (2004). 



